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1.  Introduction 


In  this  paper  we  present  a  very  fast  method  for  generating 
exponential  random  variables  in  a  digital  computer.  The  method  is 
exact,  in  the  sense  that  in  theory  it  returns  a  random  variable  with 
exactly  the  exponential  distribution.  In  practice  the  result  is  an 
approximation,  but  the  accuracy  of  the  approximation  depends  only  on 
the  word  length  of  the  computer.  The  method  is  based  on  the  general 
techniques  described  in  [1]  and  [2]  and  the  special  method  for  an 
exponential  distribution  given  in  [3] . 

When  programmed  for  the  IBM  7090  computer,  the  procedure  takes 
86  microseconds  to  generate  one  exponential  random  variable  and  requires 
390  storage  locations.  On  the  IBM  1620  computer  the  time  to  generate 
one  random  variable  is  10.8  milliseconds  and  1600  storage  locations  are 
used.  This  computer  has  a  variable  word  length,  and  a  ten  digit  number 
occupies  10  storage  locations.  The  figures  given  above  are  for  a  program 
having  ten  digit  accuracy.  We  might  remark  that  multiplication  of  two 
ten  digit  numbers  requires  17.8  milliseconds  on  this  computer. 

The  times  quoted  above  are  for  standard  FORTRAN  function  subprograms, 
and  include  linkages,  setting  index  registers,  returning  x  in 
normalized  floating  point,  etc. 

2.  Outline  of  the  procedure. 

The  procedures  for  decimal  and  binary  machines  are  essentially  the 
same.  Only  the  various  constants  differ.  The  constants  for  both  cases 
are  listed  in  Section  4. 
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We  assume  that  we  have  available  a  sequence  ui*u2****  of 
independent  uniform  [0,1]  random  variables.  The  problem  is  to 
generate  an  exponential  random  variable  x  in  terms  of  the 


Let  f  be  the  exponential  density 
f(t)  =  e_t,  0  <  t. 

Just  as  in  [2]  we  write  f  as  a  mixture  of  three  densities! 
f(t)  =  r1g1(t)  +  r2g2(t)  +  r^g^Ct). 
rigi(t)  =  e  for  0  <  kc  <  t  <  (k  +  l)c  <  4. 

r  g  (t)  =  e  ^  -  e  for  0  <  kc  <  t  <  (k  +  l)c  <  4. 

d.  (1  ~  — 

r^g^(t)  =  e  t  for  4  <  t. 

Here  k  takes  only  integral  values  and  the  constant  c  is  .1  for 
a  decimal  machine  and  .0625  for  a  binary  machine. 

Our  procedure  for  generating  x  is  essentially  as  follows. 
Generate  a  uniform  [0,1]  random  variable  u.  If  u  <  r^,  generate  a 
random  variable  with  density  g^.  If  r^  <  u  <  r^  +  r 2,  generate  a 
random  variable  with  density  g2#  If  r^  +  r2  <  u  ,  generate  a  random 
variable  with  density  g^. 
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Examination  of  the  flow  chart  will  reveal,  however,  that  the 
tests  on  u  do  not  occur  in  exactly  the  above  order,  and  also  that 
they  are  combined  with  other  tests  used  to  generate  the  random 
variables  with  densities  g^  and  g^. 

A  random  variable  with  density  g^  may  be  generated  as  y^  +  cu, 
where  y^  is  discrete,  and  u  is  uniform  on  [0,1],  The 
discrete  random  variable  y^  assumes  values  0,c,2c,....  It  is 
generated  by  the  method  described  in  [1], 

A  random  variable  with  density  g^  is  produced  as 
y^  +  c  min( v^ , v^ , . . . ,v^) ,  where  y^  is  discrete,  taking  values 
0,c,2c,....  The  v^  are  independent  uniform  [0,1],  and  z  is  a 
discrete  random  variable  with  the  distribution 

Prob(z  =  k)  =  ck/[k!(ec  -  1  -  c>] 
for  k  =  2,3,.... 


Finally  we  can  generate  a  random  variable  with  density 
sum  w  +  s,  where  s  has  density  ( r^g^  +  r2®2^/^ rl  +  T2^  ’ 
is  discrete  with  the  distribution: 


gj  as  a 
and  w 


-hk  4 

Prob(w  =  4k)  =  e  (e  -  1),  k  =  1,2,.... 
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The  random  variable  s  is  generated  by  using  the  procedure 
to  generate  the  exponential  random  number  x  with  a  slight 
modification.  To  be  exact  we  generate  an  independent  uniform 
[0,1]  random  variable  v  and  set  u  =  (r^  +  r2)v.  The 
procedure  to  generate  x  is  re-entered  at  the  first  test  on  u, 
but  now  using  u  =  (r^  +  r2)v,  which  is  uniform  on  [0,r^  +  r2^* 
This  results  in  the  generation  of  a  random  variable  with  density 

^rlSl  +  r2S2^rl  +  r2^* 

3.  FLow  chart 

In  this  flow  chart  "generate  u"  means  generate  a  uniform 
random  variable  u.  The  number  r  is  assumed  to  have  the 
fractional  representation  .d^d2.,.,  where  the  representation  is  in 
octal  if  the  program  is  for  a  binary  space  machine.  With  this 
notation  d^d2  is  nn  integer,  the  representation  being  in  octal  for 
a  binary  machine.  The  number  6  is  1/10  for  a  decimal  machine  and 
1/16  for  a  binary  machine. 
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4.  Constants 


For  both  binary  and  decimal  machines  and  are  given  by 

the  following  formulas.  The  constant  c  is  .1  for  a  decimal  machine 


and  .0625  for  a  binary  machine. 


a,  pk.i- 

e-4(k+l) 

k  =  0,1,.. 

• 

(2)  Q,  =  E  . 
k  j=2 

=3/Cj:(ec  -  1  -  c)] 

k  =  2,3,.. 

• 

The  S.  and 

X 

n.  are 

1 

Decimal 

Binary3' 

S1 

.77 

.458 

S2 

.916 

•7108 

S3 

.9314 

.7427g 

nl 

198 

294 

n2 

-764 

-144 

n3 

-9074 

-3639 

The  B.  and 

1 

C.  are 

1 

Decimal 

Machine 

Binary 

Machine 

i 

B 

C 

B 

C 

0 

.9361625820 

.931^837418 

.94510842161 

.94323480082 

1 

.9393519683 

.9362156480 

.94628827215 

.94528539525 

2 

.9415112727 

.939^089343 

.94745642249 

.94651430670 

3 

.9430881445 

.9415306939 

.94826073606 

.94757146958 

4 

.9442217866 

.9431777963 

.9489^248899 

.94843821076 

5 

.9447661013 

.9442476710 

.94932953931 

.94907597231 

6 

.9490727662 

.9448391766 

.94958329340 

.94951504666 

7 

.9529852986 

.9491545883 

.94980003626 

.94967253471 

8 

.9564834733 

.9530173032 

.94992008223 

.94985597045 

9 

.9594024969 

.9565646369 

.95190460176 

.95014449802 

10 

.9620288292 

.9594610273 

.95361457369 

.95196110918 

11 

.9643851689 

.9620617256 

.95525892439 

.953/0563829 

12 

.9663753874 

.9644731130 

.95686066511 

.95548989495 

13 

.9681837449 

.9664624958 

.95839152240 

.95693234505 

14 

.9696461638 

.9682369242 

.95974907174 

.95846135215 

The  subscript  8  means  the  number  is  given  in  octal  representation. 
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Decimal  Machine  Binary  Machine 


i 

B 

C 

B 

C 

15 

.9709809929 

.9697058602 

.96109651611 

.95996010769 

16 

.9721477967 

.9709940089 

.96237258589 

.96116288541 

17 

.9731607906 

.9722161491 

.96344671642 

.96237915952 

l8 

.9740454265 

.9731906794 

.96447806840 

.96359303267 

19 

.9748756955 

.9741022884 

.96548445 255 

.96465303841 

20 

.9756090315 

.9749092238 

.96632678306 

.96554574172 

21 

.9762878864 

.9756546743 

.96710374986 

.96637002937 

22 

.9769411563 

.9763682022 

.96786063548 

.96732383419 

23 

.9774820467 

.9770129516 

.96853524877 

.96788774279 

24 

.9779150014 

.9774905466 

.96920558798 

.96873186249 

25 

.9783264218 

.9779423592 

.96986129425 

.96928987221 

26 

.9787218719 

.9784074281 

.97048802417 

.96987974853 

27 

.9790899374 

.9787424232 

.97106734407 

.97064928279 

28 

.9794260891 

.9791686442 

.97158600948 

.97126042295 

29 

.9797130875 

.9795023095 

.97208825475 

.97180092567 

30 

.9799999300 

.9797154095 

.97258093545 

.97218820324 

31 

.9802789666 

.9800882467 

.97305557876 

.97261055489 

32 

.9805168326 

.9802838868 

.97352084773 

.97317426279 

33 

.9807267300 

.9805541596 

.97397651360 

.97360757587 

34 

.9809269259 

.9807990826 

.97442652331 

.97421630891 

35 

.9811028125 

.9809466642 

.97485411795 

.97458419727 

36 

.9812764735 

.9811351847 

.97525425825 

.97503048612 

37 

.9814292281 

.9813135507 

.97563794313 

.97545242949 

38 

.9815580886 

.9814534192 

.97601648897 

.97577828494 

39 

.9816843611 

.9815896525 

.97638555010 

.97607969003 

40 

.97672015684 

.97656636068 

41 

.97705105599 

.97695777390 

42 

.97735451196 

.97727717831 

43 

.97763681248 

.97743933438 

44 

.97791811497 

.97781881678 

45 

.97819899430 

.97816010843 

46 

.97847642598 

.97833194781 

47 

.97874079516 

.97856652131 

48 

.97899467055 

.97894151997 

49 

.97923892352 

.97913322108 

50 

.97947750310 

.97941727563 

51 

.97971601551 

.97958029084 

52 

.97993661346 

.97985429235 

53 

.98015717312 

.98004465340 

54 

.98036825703 

.98024848045 

55 

.98057295172 

.98053642199 

56 

.98077711684 

.98073305340 

57 

.98094419133 

.98078047612 

58 

.98109613615 

.98102348810 

59 

.98124587954 

.98119594925 

60 

.98139402142 

.98130639106 

6l 

.98151460394 

.98145802551 

62 

.98163244058 

.98159104690 

63 

.98168436m 

.98163745598 
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For  a  decimal  machine  the  A. 


i 

A. 

i 

A. 

i 

1 

1 

0 

0.0 

1 

.4 

2 

5 

2.2 

6 

.1 

7 

10 

.6 

11 

.7 

12 

15 

1.3 

16 

1.4 

17 

20 

1.9 

21 

2.0 

22 

25 

2.5 

26 

2.7 

27 

30 

2.8 

31 

3.2 

32 

35 

3.4 

36 

3.5 

37 

40 

.2 

4l 

.2 

42 

45 

.3 

46 

.3 

47 

50 

.6 

51 

1.2 

52 

55 

1.6 

56 

1.6 

57 

60 

2.3 

61 

2.3 

62 

65 

2.3 

66 

2.3 

67 

70 

2.4 

71 

2.4 

72 

75 

2.7 

76 

2.7 

77 

80 

3.1 

81 

3.1 

82 

85 

3.8 

86 

.5 

87 

90 

.6 

91 

.6 

92 

95 

.7 

96 

.7 

97 

100 

.9 

101 

.9 

102 

105 

1.2 

106 

1.2 

107 

110 

1.4 

111 

1.6 

112 

115 

1.7 

116 

1.7 

117 

120 

1.8 

121 

1.9 

122 

125 

2.4 

126 

2.5 

127 

130 

2.6 

131 

2.6 

132 

135 

2.9 

136 

2.9 

137 

140 

2.8 

141 

2.8 

142 

145 

3.0 

146 

3.0 

147 

150 

3.5 

151 

3.7 

152 

155 

.8 

156 

1.3 

157 

160 

1.8 

161 

1.8 

162 

165 

2.6 

166 

2.9 

167 

170 

2.9 

.  171 

2.9 

172 

175 

3.2 

176 

3.2 

177 

180 

3.6 

181 

3.6 

182 

185 

3.5 

186 

3.5 

187 

190 

3.9 

191 

3.9 

192 

195 

3.9 

196 

3.9 

197 

200 

.0 

201 

.0 

202 

205 

.1 

206 

.1 

207 

210 

.4 

211 

.4 

212 

215 

■  .4 

216 

.5 

217 

220 

.5 

221 

.7 

222 

225 

.8 

226 

.8 

227 

230 

.9 

231 

1.1 

232 

235 

1.4 

236 

1.5 

237 

240 

.0 

241 

.0 

24  2 

245 

.2 

246 

.2 

247 

250 

.2 

251 

.2 

252 

255 

.3 

256 

.3 

257 

260 

.6 

261 

.6 

262 

265 

1.0 

266 

1.1 

267 

270 

1.5 

271 

1.6 

272 

are : 


A. 

1 

i 

A1 

i 

A. 

1 

.8 

3 

1.1 

4 

1.5 

.2 

8 

.3 

9 

.5 

.9 

13 

1.0 

14 

1.2 

1.6 

18 

1.7 

19 

1.8 

2.1 

23 

2.3 

24 

2.4 

2.6 

28 

2.9 

«29 

3.1 

3.0 

33 

3.3 

34 

3.6 

3.7 

38 

3.8 

39 

3.9 

.2 

43 

.3 

44 

.3 

.3 

48 

.3 

49 

.6 

1.2 

53 

1.2 

54 

1.2 

1.6 

58 

1.6 

59 

1.6 

2.3 

63 

2.3 

64 

2.3 

2.3 

68 

2.4 

69 

2.4 

2.4 

73 

2.5 

74 

2.5 

2.7 

78 

2.7 

79 

2.7 

3.1 

83 

3.4 

84 

3.4 

.5 

88 

.5 

89 

.6 

.6 

93 

.6 

94 

.6 

.7 

98 

.9 

99 

.9 

.9 

103 

1.0 

104 

1.0 

1.3 

108 

1.3 

109 

1.3 

1.6 

113 

1.7 

114 

1.7 

1.7 

118 

1.8 

119 

1.8 

1.9 

123 

2.0 

124 

2.4 

2.5 

128 

2.5 

129 

2.5 

2.6 

133 

2.6 

134 

2.6 

2.9 

138 

2.8 

139 

2.8 

3.2 

143 

3.2 

144 

3.0 

3.3 

148 

3.3 

149 

3.6 

.7 

153 

.7 

154 

.8 

1.8 

158 

1.8 

159 

1.8 

1.9 

163 

1.9 

164 

2.6 

2.9 

168 

2.9 

169 

2.9 

2.8 

173 

3.2 

174 

3.2 

3.0 

178 

3.0 

179 

3.3 

3.6 

183 

3.5 

184 

3.5 

3.5 

188 

3.5 

189 

3.7 

3.9 

193 

3.9 

194 

3.9 

3.9 

198 

.0 

199 

.0 

.1 

203 

.1 

204 

.1 

.1 

208 

.1 

209 

.1 

.4 

213 

.4 

214 

.4 

.5 

218 

.5 

219 

.5 

.7 

223 

.7 

224 

.7 

.8 

228 

.8 

229 

.9 

1.3 

233 

1.3 

234 

1.4 

1.8 

238 

1.9 

239 

2.0 

.0 

243 

.0 

244 

.0 

.2 

248 

.2 

249 

.2 

.3 

253 

.3 

254 

.3 

.3 

258 

.6 

259 

.6 

.9 

263 

1.0 

264 

1.0 

1.1 

268 

1.2 

269 

1.2 

1.7 

273 

2.1 

274 

2.2 
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For  a  binary  machine  the  are: 


0 

0. 

1 

0.6250 

2 

0.6875 

3 

1.0000 

4 

1.3125 

5 

6 

3.3125 

7 

2.6875 

8 

3.3750 

9 

0.0625 

10 

0.1250 

11 

12 

0.3125 

13 

0.2500 

14 

0.3750 

15 

0.5000 

16 

0.4375 

17 

18 

0.7500 

19 

0.8125 

20 

0.8750 

21 

0.9375 

22 

1.2500 

23 

24 

1.3750 

25 

1.1875 

26 

1.1250 

27 

1.5000 

28 

1.7500 

29 

30 

1.5625 

31 

1.4375 

32 

1.6875 

33 

1.6250 

34 

2.1875 

35 

36 

2.1250 

37 

2.3125 

38 

2.0625 

39 

1.8125 

40 

2.5000 

41 

42 

3.1875 

43 

2.2500 

44 

2.9375 

45 

3.8750 

46 

2.5625 

47 

48 

3.5625 

49 

2.8750 

50 

3.4375 

51 

2.6250 

52 

3.1250 

53 

54 

2.7500 

55 

3.9375 

56 

3.7500 

57 

2.4375 

58 

3.2500 

59 

60 

3.0625 

61 

3.5000 

62 

3.8125 

63 

3.6875 

64 

0.1250 

65 

66 

0.1875 

67 

0.1875 

68 

0.1875 

69 

0.1875 

70 

0.1875 

71 

72 

0.9375 

73 

0.9375 

74 

1.2500 

75 

1.2500 

76 

1.2500 

77 

78 

1.1250 

79 

1.1250 

80 

1.1250 

81 

1.6250 

82 

2.3125 

83 

84 

2.3125 

85 

2.3125 

86 

2.0625 

87 

2.0625 

88 

1.8125 

89 

90 

3.0000 

91 

2.9375 

92 

2.9375 

93 

3.8750 

94 

3.8750 

95 

96 

2.3750 

97 

2.3750 

98 

2.3750 

99 

3.5625 

100 

3.5625 

101 

102 

3.5625 

103 

3.5625 

104 

2.8750 

105 

2.8750 

106 

2.8750 

107 

108 

3.4375 

109 

3.4375 

110 

3.4375 

111 

3.4375 
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